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1. INTRODUCTION 



^ I Correlated Basis Function (CBF) theory is a powerful tool to investigate the 

properties of many-body strongly interacting systems in several fields of physics. Its 
first succesful applications are found in homogeneous systems such as liquid Helium, 
electron systems (both in the forms of electron fluids and lattice structures) and 
^ ' infinite nuclear and neutron matter. The flexibility of CBF based approaches results 

Ch ■ in a realistic description of ground state (energy, momentum distribution, distribution 

functions and so on) as well as of dynamical (cross sections) quantities. While the 
I g.s. is often better studied by MonteCarlo (MC) methods, CBF and perturbative 

' theories in a CBF basis still provide one of the very few viable ways for an affordable, 

quantitative study of the response of many-body systems to external probes. In 
addition , the nuclear interaction and the induced nucleon-nucleon correlations are 
strongly dependent on the states of the nucleons themselves, and the MC techniques 
adopted, for instance, in the case of liquid Helium g.s. are not easily extended to 
nuclear cases. 

As most of the world we know is not homogeneous, CBF addicted began in the 
last decade to extend the theory (and the connected technologies) to such interest- 
ing objects as Helium films and droplets and finite nuclear systems (usually named 
nuclei). In Helium the interaction is known and has a relatively simple form, so 
the difficulty lies in the large densities. On the contrary, nuclei are considerably 
less dense and the actual problem is the aforementioned complicated structure of 
the interaction. Presently, there is a number of methods that exactly solve the light 
nuclei Schrodinger equation for realistic hamiltonians. Faddeev [1,2], Green Function 
Monte Carlo (GFMC) [3] and Correlated Hyperspherical Harmonics Expansion [4] 
are used for the A=3,4 nuclei, and GFMC has been recently pushed up to A=7 [5]. 
Variational MonteCarlo (VMC) [6] methods has been also used in light nuclei and. 
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if the spanned variational wave function space is large enough, then its description 
may be quite accurate (even if not exact). VMC has been extended to heavier nuclei, 
as 1*^0 [7]. 

As the mass number. A, grows larger and the region of medium-heavy nuclei 
is approached, using MC with realistic hamiltonians becomes increasingly difficult 
and CBF theory and cluster expansions represent an alternative and a competitive 
approach. In a series of papers, doubly closed shell nuclei (both in Is and jj coupling) 
were studied using the Fermi hypcrncttcd chain (FHNC) [8] summation technique 
[9,10,11]. Semirealistic, central interactions and simple two-body correlations, de- 
pending on the interparticle distances and, at most, on the isospin of the correlated 
pair, were used and nuclei ranging from ^He to ^*^^Pb were investigated. In a recent 
paper the same method has been extended to interactions and correlations containing 
spin, isospin and tensor components for ^^O and '^°Ca nuclei, having doubly closed 
shells in Is coupling [12]. The review of the results obtained in this work will be the 
object of the present contribution. 

We consider a standard nuclear hamiltonian 

containing only a two-nucleon potential, Vjj, whose large interparticle distances be- 
havior is dominated by meson exchange processes (one pion exchange, OPE). The 
intermediate and short distances parts of the potential are generally treated in semi- 
microscopic or purely phenomenological ways and several recipies are on the market. 
Their common feature is the accurate fit to the large body of available two nucleon 
scattering data up to energy ~ 350 MeV. For our purposes, we shall use a truncated 
version of the realistic Urbana fi4 model (U14) [13]. U14 is pametrized as the sum 
of 14 components, 

V1A,^J = E ^'(^^.OOf. ' (2) 
p=l,14 

with 

01=^^^ = [l,ai-aj,Si^]®[l,n-T^] , (3) 

Sij is the tensor operator and the remaining p > Q components contain linear and 
quadratic spin-orbit (L • S) and terms. The p > 6 components are not retained 
in the vq (U6) truncation. 

A CBF g.s. correlated A-nucleon wave function can be written as 

^{l,2...A) = [S\[F,,,]^l,2...A) , (4) 

i<i 

i.e. a symmetrized product of two-body correlation operators, F^-, acting on a mean 
field state, $(1, 2. ..A), given by a shell model wave function built up with (j)o.{i) single 
particle wave functions. is chosen of a form consistent with the interaction, 

Fi3 = E fi^^^Mj ■ (5) 

p=l,6 
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The correlation functions /^(r) are variational since they contain a set of param- 
eters fixed by minimizing the g.s. expectation value of the hamiltonian, (H) = 
(\E'|i/ |\E')/(\E'|\E'). (H) is expanded in Mayer- like cluster diagrams and infinite classes 
of these diagrams are summed by FHNC integral equations. The cluster expansion 
method has been widely applied to both finite and infinite, interacting systems with 
state independent (Jastrow) correlations . 

For a realistic approach to nuclear systems, the correlation operators are strongly 
state dependent and do not commute between each other. This fact prevents from 
the development of a complete FHNC theory for the correlated wave function of 
eq.(4). A single operator chain (SOC) approximation was adopted in ref.[14] for 
the operatorial [p > 1) correlations, together with a a full FHNC treatment of the 
Jastrow {p = 1) part in nuclear matter (NM). FHNC/SOC apparently provides a 
satisfying description of infinite nuclear matter at saturation density and a neutron 
matter equation of state compatible with the neutron stars observational data [15]. 
However, no exact check for FHNC/SOC is presently available in infinite nucleon 
matter, apart the evaluation of some additional classes of diagrams. The estimated 
accuracy in the g.s. energy has been set to less than 1 MeV at saturation [16,15]. 

We shaU use FHNC/SOC theory to study the g.s. of ^^O and 4°Ca with the U6 
interaction and the ^^O results will be compared with cluster MC (CMC) results. In 
the CMC method the Jastrow contribution is exactly treated by MC sampling, and 
that from the operatorial components is approximated by considering (via MC) up 
to four-body cluster terms. Higher order contributions are then extrapolated. 

2. ONE- AND TWO-BODY DENSITIES AND 

THE EXPECTATION VALUE OF THE HAMILTONIAN 

The one- and two-body densities, pi(ri) (OBD) and p2(i"i)i'2) (TBD), enter 
(more or less explicitely) the evaluation of (H). They are defined as 

Pi(ri) = (E^(ri-r»)) , (6) 

i 

and 

Pl(ri,r2) = (5^5(ri-r,)5(r2-r,)Of^.) . (7) 

The FHNC theory for the TBD in presence of Jastrow correlations is described at 
length in ref.[9] and here it is given as granted. Let's just recall that it is written 
in terms of the Jastrow correlation, /"^(r) = f^{r), and of the nodal (or chain) and 
elementary (or bridge) functions, Nj:y{ri, r2) and Ej:y{ri, r2), representing the sums of 
the diagrams having those topological structures, respectively. The (xy) classification 
corresponds to the exchange character of the external points (1,2): x{y) = d,e with 
d ^direct (the point does not belong to any exchange loop), e =exchange (the point 
does belong to a closed exchange loop). Futhermore, [xy = cc) (c =cyclic) diagrams 
are present, whose external points both belong to the same, non closed exchange 
loop. 
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In presence of operatorial correlations, the nodal functions acquire a state depen- 
dence, NPy{ri,r2)- Because of the non commutativity of the correlations, the same 
topological diagrams may originate from several different orderings. Spin-isospin 
traces provide a weight to each of the diagrams, which, in turn, depends on the or- 
dering itself. Keeping track of them is not feasible by FHNC-like integral equations 
and only selected classes of diagrams may be correctly summed. The SOC approx- 
imation consists in summing p > 1 chains, where each link may contain just one 
operatorial element and Jastrow dressings at all orders. We remind that operatorial 
dependence comes also in account of the exchange of two nucleons, as the exchange 
operator is i^^ = -[1 + Oj + Oj + Ofj]/4. 

The FHNC/SOC equations are given in ref.[12], where they are solved in the 
FHNC/0 approximation (corresponding to set to zero the elementary diagrams). 
The validity of this seemingly crude approximation will be discussed later. 

The structure of the OBD is 

Pi(ri) = pf(ri) [1 + C/7(ri)] + C/°"(ri)Cd(ri) , (8) 

with 

pf(ri)= [po(ri) + t//(ri)]Cd(ri) . (9) 

Po(ri) = Eal<^a(l)l^ is the mean field density, Cd(ri) = exp {C//(ri)}, t/'/(g)(ri) 

are Jastrow vertex corrections and U^^^-^i'^i) represent operatorial vertex corrections, 

linked, by integral equations given in refs.[9,12], to the FHNC nodal functions. 

A measure of the accuracy of the computed OBD is the degree of fulfilment of 
its normalization. Si = A, where 

Si= J dVipi(ri) . (10) 

In addition, the exact TBD has to comply with the pair saturation, S2 — I, and (for 
the nuclei we are considering) the isospin saturation, Sj- = —1, properties, with 

^2 = j^^^ I d'r, I dV2pr'(ri,r2) , (11) 

and 

Sr = ^J d'n J d\2Pr\r,,r2) . (12) 

A spin saturation sum rule, = —1, holds only in absence of tensor correlations. 

Deviations of the sum rules from their exact values are due to (z) the FHNC/0 
scheme and (ii) the SOC approximation. In ref.[9] was found that E^^, i.e. the sum 
of the ee-elementary diagrams whose external points belong to the same exchange 
loop, may substantially contribute to both and to the potential energy, if the 
potential has large exchange terms. 

In evaluating (H), it is convenient to use the Jackson-Feenberg identity [17] for 
the kinetic energy, (T) , with the result 

{T)^Tjp^T^ + Tf , (13) 
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Figure 1. Correlation functions for Ca and NM at saturation density. The oper- 
atorial correlations (right panel) are given as lines for ^°Ca and symbols for NM. 



with 

T^ = -A—{^*G^Vl^-{Vi^*)G^{Vi^)) , (14) 

and 

T^ = -A— ($* [GVfG'- ViC- ViGl $) , (15) 
4m 

where G = S]\F,j. Ref.[12] describes how to compute in FHNC/SOC. 

To evaluate Tp and the two-body potential energy {v) = V2, we define an inter- 
action energy W — Tp + V2 = {Hjp)- The two-body operator, H^jp{ri2)i is 

H'j^{r^2) = -^5,1 {r(ri2)VV'(ri2) - ^^fin^) ■ V/'=(ri2)} 

+firi2)v^ir,2)fHri2) • (16) 

In FHNC/SOC, W is split into four parts, W = Wq + Ws+Wc + W^s, where Wq is 
the sum of the diagrams with only p = 1, central chains between the interacting points 
(IP), connected by Hjp; Wg sums diagrams having operatorial vertex corrections 
reaching the IP's and central chains; Wc contains diagrams with one SOC between 
the IP's; Wcs contains both. The lenghty formulas for W are again found in ref.[12]. 



3. ENERGETICS 

The results presented in this section have been obtained with the single particle 
wave functions, </>q:(^)) generated by a harmonic oscillator well with parameter b. 
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Table 1 



Breakup of the energy per nucleon, in MeV, for O with the U6 potential. The J 
column gives the energies with the Jastrow corelated wave function, the O column 
those with the operatorial correlation and the CMC columns the cluster MonteCarlo 
results of ref.[18]. Egs = {{H) - Tcm)/A. 





J 


CMC(J) 





CMC(O) 




0.88 


0.93(28) 


2.33 


2.35(43) 




1.25 


1.27(08) 


1.97 


2.00(13) 








0.25 


0.27(01) 




2.40 


2.43(12) 


2.29 


2.23(14) 


iv') 


-26.59 


-26.24(26) 


-32.03 


-30.12(42) 








-10.28 


-9.77(09) 


{v) 


-22.07 


-21.56(25) 


-35.47 


-33.03(31) 


(T) 


24.61 


24.33(21) 


31.16 


29.45(33) 


(H) 


2.54 


2.77(09) 


-4.33 


-4.59(10) 


Egs 


1.72 




-5.15 





We kept b fixed (6(^^(9) = 1.543 fm and 6(4°Ca) = 1.654 fm), even though it 
might be considered as a variational parameter. As far as the variational correlation, 
Fij, is concerned, we follow ref.[9] and adopt an Euler correlation operator obtained 
by minimizing the energy at the second order of the cluster expansion, {H2)- A 
corresponding choice was done in ref.[15] for NM. The variational parameters of the 
functions /^(r) arc the /lea/m^' distances, dp, where /^(r > di) = 1, fP^^{r > dp) = 
and {f^)'{r = dp) = 0. As in NM, we use only two healing distances: dc and dt, for 
the central {p — 1, 2, 4, 5) and tensor {p — 3, 6) channels, respectively. An additional 
parameter is the quenching factor, a, of the p > 1 components of potential in the 
Euler equations. We have taken the U14 NM parameters, dc = 2.15 fm, dt = 3.43 fm 
and a = 0.8, on the assumption that the short range behavior of the correlation does 
not vary drastically in going from the infinite to the finite systems. It is clear that a 
more accurate variational search for each nucleus would provide a lower ground state 
energy. 

The ^°Ca correlation functions are shown in Figure 1 and compared with the 
corresponding NM functions. They are similar, especially the longer ranged tensor 
parts. The most visible differences are found in the spin and isospin components and 
in the shortest range part of 

The error in the Si sum rule, in FHNC/0, is well less than 1% in both nuclei 
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for either the Jastrow (J) or the operatorial (O) correlations. The same accuracy is 
found for S2 and Sj- in the J-model after the insertion of the E^^ diagram (FHNC-1). 
The situation worsens for the 0-model, where FHNC/SOC violates the sum rules by 
a maximum amount of ~ 9%, as it was already found in NM [15]. 

The U6 energetics is displayed in Table 1 for -^^O and Table 2 for ^°Ca, for the 
J and O correlations. The Tables provide the potential, kinetic and total energy per 
nucleon. The separated expectation values of the components of the potential are 
also shown, together with the ground state energy. Egg, obtained by subtracting the 
center of mass kinetic energy, Tcm, from (H). 

The ^^O results are compared with the CMC calculations of Pieper [18]. FHNC 
shows an error of 1 — 2% for (T) and (^2) in the J case. The total energy error is 
larger (~ 9%) as (H) is given by the cancellation of two large numbers. The same 
situation is met in the 0-model: the kinetic and potential energy errors are 5 — 7%. 
The absolute error in {H) is well less than 1 MeV, consistently with the estimated 
accuracy of FHNC/SOC in NM [15]. The OPE parts of the potential provide most 
of the binding since, in its absence, ^^O is not bound. The same holds in ^^Ca, where 
the introduction of tensor correlations and potentials increases the kinetic energy by 
~ 5.6 MeV, compensated by an additional potential energy contribution of ~ —13.6 
MeV, providing a bound nucleus. 

In order to compare the theoretical energies with the experimental values one has 
to estimate the mean value of the full U14 interaction. Its p > 6 components represent 
the momentum dependent (MD) part of the potential, whose expectation value is 
more difficult to evaluate than U6 by means of cluster expansions. The usual strategy 
in NM is to include only low order (second and third) cluster MD contributions to 
(v^^). This procedure is probably accurate for U14, since its MD part contributes 
to the g.s. energy by less than 1 MeV at piVM = 0.16 fm~^ i^Ef^M = 0-44 McV), 
but it could be questionable for other potentials, as the recent Argonne vis model 
[19]. Moreover, it is known that a pure two-body interaction does not reproduce the 
NM saturation (energy, density and incompressibility) and the binding energies of 
light nuclei, so three- nucleon interactions, Vijk, are to be introduced. The TNI model 
of ref.[13] approximated the effect of Vijk by adding two density dependent terms to 
U14: a repulsive TNR one, reducing the effect of the intermediate range part of Vij, 
and an attractive TNA part. The TNR term was taken as 

UU + TNR^ Yl + ^/(^^^O + ^sir^JmJ , (17) 

p=l,14 

where j g are the long-, intermediate- and short-range parts of the potential, Pnm 
is the NM density and 71 = 0.15 fm~^. The TNA term was assumed to contribute 
as 

= 372p?,Me"^^''^^ , (18) 
with 72 = —700 Mev fm^ and 73 ~ 13.6 fm"^. The values of the 7-parameters were 
obtained by fitting the NM saturation properties in FHNC/SOC. 

Here we evaluate the TNR contribution by using [pi(ri)pi(rj)]^/^ in place of 
Pnm in (17) and 

TNA = ^ J d^npiin) 3-f2pi{ri)e-^'P'^'''^ . (19) 
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Table 2 

As Table 1 for ^°Ca. The CMC results are not available in this case. 





J 







-1.41 


-0.21 




1.57 


2.30 






0.29 




2.99 


2.71 


iv') 


-32.40 


-39.48 






-14.14 


{v) 


-29.26 


-47.28 


(T) 


30.55 


39.69 


(H) 


1.30 


-7.59 


Egs 


1.01 


-7.87 



Table 3 

Energies per nucleon (in MeV) for ^^O and ^°Ca with the U14-I-TNI interaction. See 
text. 





{H)u6 




TNA 


{H)ui4+TNI 


AE, 


Egs 


fpexpt 


160 


-1.88 


0.33 


-4.53 


-6.08 


0.05 


-6.85 


-7.72 


40Ca 


-3.20 


0.68 


-4.59 


-7.11 


0.05 


-7.35 


-8.30 



AE is estimated in local density approximation (LDA) as 

A£;^^ = I / rf'npi(ri)A<^[pi(ri)] . (20) 

Table 3 gives the ^^O and ^'^Ca energies with the U144-TNI interaction and the 
operatorial correlation model. AE, is the Coulomb energy correction taken from 
ref.[ll]. The experimental ground state binding energies per nucleon are underesti- 
mated by ~ 1 MeV in both nuclei. The MD contribution is enough low to makes us 
enduring the LDA in its evaluation. We are considerably satisfied with these results, 
most af all in consideration of the further energy lowering that could be obtained by 
a careful variational minimization and of the FHNC/SOC accuracy. 
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4. CONCLUSIONS 



In this contribution we have presented some resuhs for the ground state energy 
of doubly closed shell nuclei in Is coupling (^^O and ^'^Ca) in the framework of the 
correlated basis function theory. The correlated wave function includes central and 
tensor correlations, whose non commutativity does not allow for a complete FHNC 
treatement, contrary to the purely scalar case. We have extended the single opera- 
tor chain approximation scheme (as used in realistic studies of nuclear and neutron 
matter) to the finite nuclei. The U6 truncation of the realistic Urbana V14 nucleon- 
nucleon potential has been adopted as a test. 

By the analysis of the sum rules we have shown that FHNC/SOC provides 
an accurate one-body density. A comparably good accuracy is obtained for the 
normalizations of the two-body density (^'2 and Sr) when tensor correlations are not 
included. Their insertion slightly worsens the excellent fuUfiment of S2,t, but the 
maximum violation is ~ 9 %, similar to what was found in nuclear matter. 

The comparison of the energy contributions in ^^O for the U6 interaction with 
the cluster Monte Carlo results shows a maximum disagreement varying from ~ 2 % 
for the Jastrow model to ~ 7 % for the tensor model. The absolute error in the 
ground state energy per nucleon is well less 1 MeV, compatible with the estimated 
accuracy of the FHNC/SOC approach in nuclear matter at saturation density. 

The ground state expectation value of the complete realistic U14-I-TNI interac- 
tion has been evaluated by estimating the small momentum dependent contribution 
in local density approximation. The binding energies are about 1 Mev lower than 
the experimental values. Further improvements are reasonably to be expected by a 
better variational minimization. To our knowledge and opininion, this is one of the 
very few reliable microscopic evaluations of the g.s. energy of these nuclei with a 
realistic hamiltonian. Recently derived nuclear potentials provide a higher quality 
fit to the nucleon-nucleon scattering data, however their momentum dependence is 
stronger than in U14. Consequently, a better than LDA treatment of this part is 
imperative if they are to be used in the nuclei we have studied. Our group is working 
along this line, as well as to the insertion of explicit three-nucleon forces. 
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